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Abstract 

The contamination of groundwater by heavy metal ions around a lead and zinc plant has been studied. As a case 
study groundwater contamination in Bonab Industrial Estate (Zanjan-lran) for iron, cobalt, nickel, copper, zinc, 
cadmium and lead content was investigated using differential pulse polarography (DPP). Although, cobalt, copper 
and zinc were found correspondingly in 47.8%, 100.0%, and 100.0% of the samples, they did not contain these 
metals above their maximum contaminant levels (MCLs). Cadmium was detected in 65.2% of the samples and 
17.4% of them were polluted by this metal. All samples contained detectable levels of lead and iron with 8.7% and 
13.0% of the samples higher than their MCLs. Nickel was also found in 78.3% of the samples, out of which 8.7% 
were polluted. In general, the results revealed the contamination of groundwater sources in the studied zone. The 
higher health risks are related to lead, nickel, and cadmium ions. Multivariate statistical techniques were applied for 
interpreting the experimental data and giving a description for the sources. The data analysis showed correlations 
and similarities between investigated heavy metals and helps to classify these ion groups. Cluster analysis identified 
five clusters among the studied heavy metals. Cluster 1 consisted of Pb, Cu, and cluster 3 included Cd, Fe; also each 
of the elements Zn, Co and Ni was located in groups with single member. The same results were obtained by 
factor analysis. Statistical investigations revealed that anthropogenic factors and notably lead and zinc plant and 
pedo-geochemical pollution sources are influencing water quality in the studied area. 

Keywords: Heavy metals, Groundwater, Differential pulse polarography, Multivariate data analysis 



Introduction 

Water is one of essential compounds for all forms of 
plants and animals [1], thus its pollution is generally con- 
sidered more important than soil and air. Due to its spe- 
cific characteristics, this liquid bears unique properties. It 
is the most effective dissolving agent, and adsorbs or sus- 
pends many different compounds [2]. 

More than one billion people in the world do not have 
suitable drinking water, and two to three billions lack ac- 
cess to basic sanitation services. About three to five mil- 
lions die annually from water related diseases [3]. 

Surface water (fresh water lakes, rivers, streams) and 
groundwater (borehole water and well water) are the prin- 
cipal natural water resources. Nowadays one of the most 
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important environmental issues is water contamination 
[4,5]. Heavy metals are among the major pollutants of 
water sources [6]. Despite this, heavy metals are sensitive 
indicators for monitoring changes in the marine environ- 
ment. Due to human industrial activities, the levels of heavy 
metals in the aquatic environment are seriously increasing 
and have created a major global concern [7,8]. Some of 
these metals are essential for the growth, development and 
health of living organisms, whereas others are non-essential 
as they are indestructible and most of them are categorized 
as toxic species on organisms [9]. Nonetheless the toxicity 
of metals depends on their concentration levels in the en- 
vironment. With increasing concentrations in environment 
and decreasing the capacity of soils towards retaining heavy 
metals, they leach into groundwater and soil solution. Thus, 
these toxic metals can be accumulated in living tissues and 
concentrate through the food chain. 
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Cadmium is regarded as the most serious contaminant 
of the modern age [10]. Copper is classified as a priority 
pollutant because of its adverse health effects [11]. Zinc 
and iron are essential elements and are generally consid- 
ered to be non-toxic below certain levels [12]. Lead is 
not an essential trace element in any organism and has 
no known biological function. It can cause a variety of 
harmful health effects [13] and is known as a fatal neu- 
rotoxicant [14]. Excessive concentrations of cobalt can 
cause death and various compounds of nickel are car- 
cinogenic [15]. These menaces provoke the studies on 
the monitoring of these heavy metals in this chain being 
important for protection of public health. 

A variety of techniques including x-ray fluorescence 
(XRF), neutron activation analysis (NAA), inductively 
coupled plasma-atomic emission spectrometry (ICP-AES), 
atomic absorption spectrometry (AAS) and graphite fur- 
nace atomic absorption spectrometry (GFAAS) have been 
used for evaluating the heavy metal concentration in envir- 
onmental samples [16-20]. Beside their valuable charac- 
teristics, these techniques suffer from some disadvantages 
such as heavy capital cost, expensive maintenance, and in- 
sufficient sensitivity for very low concentrations of metals. 
Voltammetric methods are known as sensitive techniques 
for determination of a variety of chemical species [21]; 
among these techniques, differential pulse polarography 
(DPP) bears some advantages for accurate and precise de- 
tection and determination of trace amounts of heavy metal 
ions in environmental samples [22,23] . 

Evaluation of the contaminants resulted from excavation 
of zinc and lead mines and development of related indus- 
tries in Zanjan province-Iran and their negative environ- 
mental impacts is critical and important. Lack of a 
systematic investigation of the probable heavy metals con- 
tamination around National Iranian Lead and Zinc Com- 
pany (NILZ) in Bonab Industrial Estate (BIE), in Zanjan 
province, promotes to assess the quality of groundwater 
sources in this industrial zone. These are the main sources 
of drinking water and irrigation for a part of people who 
live around NILZ Company. In this research, DPP tech- 
nique was used to determine the concentrations of seven 
heavy metals (iron, cobalt, nickel, copper, zinc, cadmium 
and lead) in water samples and the results were compared 
with the maximum contaminant levels (MCLs) specified 
by WHO as well as Institute of Standards and Industrial 
Research of Iran (ISIRI). The multivariate statistical ana- 
lysis was conducted to categorize the metals and to distin- 
guish the source of the contaminants. 

Materials and methods 

Study area 

Zanjan province (located in north west Iran), has a large 
metalliferrous site and has been considered as a traditional 
mining region since antiquity [24]. There are still large 



reserves of lead and zinc in the area. Both mines and 
smelting units within the province present a risk of con- 
tamination of soils, plants, and surface/groundwater 
resources through dissemination of particles carrying 
metals by wind action and/or by runoff from the tailings 
[25]. Transportation of concentrated ore by trucks for 
about 110 kilometers from mines in Angouran to NILZ is 
another anthropogenic source of metal contamination, es- 
pecially along the roads. 

In this study, Bonab Industrial Estate (BIE) and its 
neighborhood was selected for detailed study. The re- 
search was focused on the environmental impacts of 
NILZ Company (36° 66/ N, 48° 48/ E) located within BIE, 
about 12 km east of Zanjan city. The NILZ Company 
was established in 1992, with a current consumption of 
about 300,000 tons of raw ore and an annual production 
of 55000 tons of Pb and Zn [26,27]. The plant is situated 
over an aquifer, which is the only source of fresh water 
available in the area, supplying a part of drinking water 
to Zanjan citizens and its neighboring areas as well as 
water used for agricultural and industrial consumptions. 
The tailings from BIE, estimated to be about 2.5 million 
tons, contain a variety of toxic elements, notably Pb, Zn, 
and Cd [26]. They are damped in the vicinity of the Es- 
tate and are exposed to wind and rain, contributing to 
soil, surface and groundwater contamination. 

Sample collection and storage 

To examine the extent of the contamination by toxic 
metals leached from tailings, 23 spring/groundwater 
samples were collected and analyzed from the studied 
area. Nineteen groundwater and four spring water sta- 
tions were selected from the site within a radius of five 
km from NILZ Company (Figure 1). 

Sampling stations were selected, taking into account 
the direction of groundwater flow (west), direction of 
prevailing winds (west and south west) and also the 
density of the population within the studied area. How- 
ever, limitations on number and distribution of sampling 
stations are set due to the spatial distribution of available 
bore wells within the studied area. Table 1 shows the lo- 
cation of sampling stations for this study. 

From each station three replicate samples were 
selected for analysis. Glassware and vessels were treated 
in 10% (v/v) nitric acid solution for 24 h and were 
washed with distilled and deionized water. The samples 
were collected in polypropylene containers, labeled and 
immediately few drops of HN0 3 (ultra pure grade) to 
pH < 2 were added to prevent loss of metals, bacterial 
and fungal growth and then stored in a refrigerator. 

Reagents and standards 

All the chemicals used in this study were mostly reagents of 
highest grades (Merck) and used without further treatment. 
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Figure 1 Location map of the studied area indicating sampling points. 



The chemicals used were: dimethylglyoxime (>99%), am- 
monia solution (25%), ammonium chloride (>99.8), acetic 
acid (>99.8), hydrochloric acid (37%), nitric acid (65%), 
pyrocatechol (>99%), and sodium hydroxide (>97%). The 
heavy metal standards were prepared from stock solutions 
of 1000 ±5 mg/L (Merck) by successive dilution with 
ultra-pure water. Polargraphic mercury was used as elec- 
trode in heavy metal determination (Merck). 

Sample digestion 

Groundwater samples were filtered through 0.45 \im fil- 
ters. To ensure the removal of organic impurities from 
the samples and thus preventing interference in analysis, 
the samples were preserved and digested with concen- 
trated nitric acid. To this end 1 mL of nitric acid was 
added to water sample in 50 mL volumetric flask. 

Sample analysis in the field 

The pH, electrical conductivity (EC) and dissolved oxy- 
gen (DO) of the samples were immediately measured at 
sampling stations by using a portable digital pH meter 
(Hach HQ 40d). Recorded pH and EC of samples varied 
in the range of 7.2 -8.3 and 326-1857 (\iS /cm) 



respectively (Table 1). The pH values of the samples 
were within the WHO range (6.5 -8.5) but those of 
ECs were below the announced value of MCL by WHO 
(1500 \iS /cm), except for samples number W6 and W8. 

Sample analysis 

Water samples were analyzed for the presence of iron, 
cobalt, nickel, copper, zinc, cadmium and lead using a 
differential pulse polarography (Metrohm 797 VA). Dis- 
solved air was removed from the solutions by degassing 
with N 2 gas (99.999%) for 5-10 min prior to each run. 
Standard addition method was used for the analysis. The 
polarography parameters are given in Table 2. Digested 
samples were analyzed in triplicate and the average con- 
centrations of metals were reported in (ig/L. 

Statistical analysis 

SPSS statistical package (Window version 18) and soft- 
ware Excel 2007 are used for data analysis. The analysis of 
the experimental data was carried out by using one-way 
ANOVA, Pearson correlation matrix, Cluster Analysis, 
Principal Component Analysis (PCA) and Factor Analysis 
(FA) methods [28,29]. Pearson correlation matrix shows a 
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Table 1 GPS location and some physical properties of sampling wells 


Site/location 


GPS location (UTM) 


pH 


EC (mS /cm) 


DO (mg /L) 


t(°C) 


Depth (m) 


Distance 1 (m) 




X 


Y 
















W1 


290458 


4054526 


7.81 


576 




6.91 


16 


60 


2418 


W2 


289862 


4052820 


7.93 


465 




7.46 


16 


75 


2958 


W3 


290843 


4051745 


7.58 


729 




8.36 


15 


50 


4349 


W4 


289758 


4052297 


7.70 


758 




7.73 


15 


35 


3308 


W5 


287326 


4051641 


8.15 


900 




7.10 


28 


spring 


3617 


W6 


290619 


4051470 


7.41 


1607 




7.51 


14 


13 


4512 


W7 


284992 


4055674 


7.86 


500 




6.40 


17 


45 


3234 


W8 


284999 


4055445 


8.23 


1857 




6.95 


18 


32 


3247 


W9 


283572 


4054398 


8.27 


826 




7.22 


27 


spring 


4736 


W10 


286973 


4056943 


7.51 


361 




8.39 


17 


150 


2203 


W11 


289883 


4054679 


7.45 


369 




8.07 


19 


150 


1777 


W12 


291848 


4052965 


7.53 


480 




8.87 


15 


150 


4352 


W13 


290659 


4059821 


7.87 


334 




9.03 


16 


spring 


5320 


W14 


291539 


4058941 


8.29 


384 




7.10 


26 


spring 


5073 


W15 


290371 


4052030 


7.52 


990 




9.36 


13 


25 


3971 


W16 


289285 


4052723 


7.45 


1157 




8.22 


14 


10 


2746 


W17 


288931 


4053010 


7.43 


1073 




7.78 


14 


150 


2344 


W18 


287475 


4054510 


7.85 


815 




9.26 


18 


20 


995 


W19 


287004 


4055727 


8.08 


564 




7.59 


14 


50 


1300 


W20 


287319 


4054582 


7.20 


326 




7.20 


15 


13 


1054 


W21 


286247 


4055985 


7.35 


850 




7.70 


15 


70 


2183 


W22 


287340 


4054826 


7.50 


671 




6.97 


15 


42 


927 


W23 


290294 


4053296 


7.80 


415 




8.90 


16 


70 


2837 



1 Distance from NILZ Company. 



probable common source of the pollutants. Cluster ana- 
lysis is used for dividing the studied metal ions into the 
similar classes with respect to their normalized concentra- 
tion level PCA is designed to transform the original vari- 
ables into new, uncorrelated variables (axes), called the 
principal components. Factor Analysis is similar to Princi- 
pal Component Analysis method except for the prepar- 
ation of the observed correlation matrix for extraction and 
the underlying theory [30]. 

The one-way ANOVA method allows testing the signifi- 
cant difference of the means. For this test each sampling 
location was selected as a group and its heavy metal con- 
centration as the corresponding variable. The ANOVA 
test requires three assumptions, i.e. the random behavior 
of the occurrence, the homogeneity of variance and the 
normal distribution behavior of the metal ions in the sam- 
ple stations. These were tested by using Runs test, Levene 
statistic and the K-S (Kolmogorov-Smirnov) methods, re- 
spectively. It is noteworthy that instead of the ANOVA 
test, one can use the Kruskal-Wallis test. The latter is a 
non-parametric test without requirements announced for 



the ANOVA test [28,29]. In this work both of the methods 
were tested for a comparison. 

The bivariate correlation procedure computes the pair 
wise associations for a set of variables and displays the 
results in a matrix. It is useful for determining the strength 
and direction of the association between two variables. The 
correlation coefficients computed by bivariate correlation 
procedure lay in the range -1 (for the cases in which a per- 
fect negative relationship exists) to +1 (for a perfect posi- 
tive relationship). A value of 0 indicates there is no linear 
relationship among the variables. For normally distributed 
variables, the Pearson method can be used to calculate the 
correlation coefficient. For normally distributed variables, 
the Pearson correlation was used for bivariate correlation, 
otherwise non-parametric Spearman method was applied. 

Cluster analysis is a method for dividing a group of 
metals into classes so that similar metals, with respect to 
variable space, are in the same class. In fact, the groups 
are not known prior to applying this mathematical ana- 
lysis and no assumption is made about the distribution 
of the variables [28,29]. 
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Table 2 Instrument operating parameters for the analysis 
of the investigated heavy metals 



Parameters Heavy metals 





Fe 1 


Co and Ni 2 


Cu, Zn, Cd and Pb 3 


Working electrode 


HMDE 


HMDE 


HMDE 


Drop size 


7 


4 


4 


Stirrer speed 


2000 rpm 


2000 rpm 


2000 rpm 


Mode 


DP 


DP 


DP 


Purge time 


300 s 


300 s 


300 s 


Deposition potential 


-300 mV 


-0.7 V 


-1.15V 


Deposition time 


60s 


90 s 


90s 


Equilibrium time 


5 s 


10 s 


10s 


Pulse amplitude 


50 mV 


50 mV 


50 mV 


Start potential 


-200 mV 


-0.8 V 


-1.15 V 


End potential 


-550 mV 


-1.25 V 


0.05 V 


Voltage step 


4 mV 


4 mV 


6 mV 


Voltage step time 


0.4 s 


0.3 s 


0.1 s 


Sweep rate 


10 mV/s 


13 mV/s 


60 mV/s 


Peak potential 


-380 mV 


-1.13, -0.97V 


-0.10, -0.98, -0.56, -0.38 V 



1 1 0 mL sample solution + 1 00 uL Catechol solution (1 M) + 1 mL phosphate 
buffer; pH =7.0, 2 10 mL. 

sample solution + 100 uL dimethylgloyoxime solution (0.1 M) + 0.5 ml NH 4 CI 
pH = 9.5, 3 10 mL. 

sample solution + 1 mL ammonium acetate buffer; pH =4.6. 



The major objective of FA is to reduce the contribution 
of less significant variables to simplify even more of the 
data structure given by PC A. This goal can be achieved 
by rotating the axis defined by PCA and constructing 
new variables, also called Varifactors [31]. PCA reduces 
the dimensionality of data by a linear combination of ori- 
ginal data to generate new latent variables which are or- 
thogonal and uncorrelated to each other [32]. The major 
objective of FA is to reduce the contribution of less sig- 
nificant variables to simplify even more of the data struc- 
ture coming from PCA. All significance statements 
reported in this study are at the P < 0.05 level. 

Results 

Extent of heavy metals contamination 

The results of analysis of target metal ions i.e. Fe, Co, Ni, Zn, 
Cd and Pb in samples from 23 studied wells are given in 
Table 3. It is noteworthy that the reported values are based 
on three replicate determinations. Table 4 is prepared in 
order to give a simple comprehensive interpretation on the 
obtained data, and to compare the concentration of the stud- 
ied metals in the samples with the MCL values reported by 
WHO and ISIRI. The results show that Fe, Co, Ni, Cu, Zn, 
Cd and Pb are detected in 100%, 47.8%, 78.3%, 100%, 100%, 
65.2% and 100% of the samples, respectively. The concentra- 
tion of metals (in Lig /L) in the samples were found in the 
range of 75.90 -339.75 for Fe; ND (not detected) -99.82 for 



Co; ND -84.15 for Ni; 6.59-65.31 for Cu; 27.79 -2227.80 for 
Zn; ND -14.87 for Cd; and 0.74 -12.45 for Pb. 

Comparison of the concentration of heavy metals 

In order to deduce the frequencies of the concentration 
of each metal in the samples, the Chi-Square test was 
applied [29]. Here, the frequency means the number of 
times a given range of concentrations occurs, and the 
Chi-Square test is used to examine whether the observed 
frequencies differ significantly from those which would 
be expected on the null hypothesis. This test indicates 
that there is no significant difference between observed 
frequencies of the heavy metals. 

The random and normal distribution assumptions were 
checked by Runs and K-S methods, respectively. Another 
requirement for applying the ANOVA test is that the var- 
iances of the groups are equivalent. Based on the statisti- 
cally verification done by Levene test, the homogeneity of 
variance was found to be significant for the samples 
(Levene statistic = 5.696, P < 0.001). Although the Levene 
statistic parameter rejects the null hypothesis, as the group 
variances are equal the ANOVA test can be yet used. Al- 
ternatively, the homogeneity and normal distribution in 
the data can be achieved by transforming the obtained 
data to another mathematically presentation which lowers 
the difference between the data. This can be achieved for 
example by using the logarithmic form of data. In 
addition, one can use a non-parametric test. This type of 
tests does not require to a homogeneity assumption. 

The ANOVA method was used under two conditions. 
In fact, although the homogeneity of the data was not 
shown, ANOVA was applied to the data. In addition, by 
transforming the data as their logarithmic form, homo- 
geneity in the observed data was achieved. The non- 
parametric method Kruskal-Wallis is based on ranks of 
the data variances. This method was used for the same 
scope as ANOVA. Both parametric and non-parametric 
methods used for comparison of the concentrations of 
heavy metals among sampling sites show a statistically 
significant difference depending on sampling locations. 

Bivariate correlations of investigated heavy metals 

To deduce the probable common source of metals in 
water samples, the bivariate correlation procedure was 
used (Table 5). This procedure computes the pair wise 
associations for a set of metals and displays the results as a 
matrix. It is useful for determining the value of association 
of the investigated metals. Because, obtained data was not 
normally distributed, Spearman method was applied. 

Classification of the investigated heavy metals by 
cluster analysis 

Cluster analysis grouped the studied heavy metals into clus- 
ters (called groups in this study) on the basis of similarities 
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Table 3 Metal contents in water samples (ng/ L) from the wells 


Sample 


Fe 


Co 


Ni 


Cu 


Zn 


Cd 




Pb 


W1 


129.30 ±10.91 


9.91 ±0.91 


84.15 ± 10.84 


44.39 ± 0.26 


196.95 ±11.40 


4.55 ± 0.25 




1 1 .78 ±2.76 


W2 


339.75 ±39.77 


60.96 ± 1 .95 


52.24 ± 4.46 


26.25 ± 2.03 


169.51 ±6.68 


14.87 ±0.86 




6.10 ±0.54 


W3 


143.99 ±20.1 6 


1 1 .07 ±1.15 


41 .59 ± 1 .46 


15.27 ±1.44 


140.53 ±13.55 


0.68 ± 0.06 




1 .26 ± 0.64 


W4 


228.71 ±4.50 


3.51 ±0.82 


42.98 ±2.41 


18.39 ±1.46 


460.1 9 ±26.21 


1 .29 ± 0.45 




1 .67 ± 0.78 


W5 


194.11 ±18.49 


1.25 ±0.15 


21.82 ± 1.61 


16.65 ±1.05 


353.92 ±59.79 


ND 




2.63 ± 0.28 


W6 


109.83 ±16.98 


ND 2 


20.09 ± 1 .23 


8.83 ± 0.27 


232.38 ±41.84 


ND 




1 .76 ± 0.42 


W7 


146.51 ±3.63 


0.43 ± 0.02 


14.12 ±0.91 


15.71 ±1.93 


582.66 ±61. 54 


3.55 ±0.37 




5.46 ±0.48 


W8 


11 6.20 ±9.88 


ND 


ND 


27.71 ±1.90 


96.54 ± 6.23 


ND 




2.51 ±0.16 


W9 


288.97 ±16.40 


1.89 ±0.31 


12.1 5 ±0.44 


65.31 ±6.00 


541.43 ±22.53 


3.41 ± 0.34 




12.45 ±0.73 


W10 


302.70 ± 20.74 


2.24 ±0.1 5 


6.1 3 ±0.22 


35.81 ±3.41 


60.36 ± 6.74 


0.36 ± 0.05 




6.52 ±0.31 


W11 


124.80 ±7.37 


ND 


5.20 ±0.21 


12.06 ±0.67 


205.36 ±12.64 


ND 




2.04 ±0.1 9 


W12 


144.21 ±15.68 


99.82 ± 9.64 


9.55 ± 0.41 


21.96 ±0.92 


84.54 ± 6.23 


0.62 ± 0.20 




5.45 ± 0.38 


W13 


132.58 ±6.49 


ND 


ND 


9.84 ±1.11 


49.26 ±2.11 


1 .23 ± 0.23 




3.25 ± 0.24 


W14 


308.73 ± 22.94 


ND 


12.21 ±1.11 


13.79 ±1.09 


43.68 ±2.1 5 


0.83 ±0.11 




6.71 ±0.24 


W15 


120.64 ±24.25 


1.52 ±0.27 


8.04 ±0.28 


60.77 ±3.10 


11 3.99 ±7.1 2 


ND 




6.56 ±0.31 


W16 


75.90 ±10.80 


3.1 9 ±0.35 


7.68 ± 0.47 


31.95 ±2.47 


90.04 ±7.71 


ND 




6.60 ± 0.24 


W17 


245.96 ±29.36 


ND 


7.21 ±0.24 


21.52 ±1.66 


11 5.69 ±15.27 


1 .43 ± 0.36 




3.31 ±0.27 


W18 


179.67 ±22.70 


ND 


6.37 ±0.27 


8.22 ± 0.72 


31.72 ±4.87 


0.73 ± 0.05 




0.91 ±0.07 


W19 


10 1.49 ±6.39 


ND 


ND 


24.87 ±0.1 2 


27.79 ±1.51 


0.59 ± 0.06 




0.74 ±0.21 


W20 


84.60 ±12.10 


ND 


ND 


28.1 7 ±1.02 


133.82 ±15.59 


0.1 5 ±0.06 




4.60 ± 0.22 


W21 


11 5.60 ±32.23 


ND 


ND 


6.59 ± 0.92 


65.99 ±4.21 


2.25 ± 0.39 




2.25 ± 0.32 


W22 


137.38 ±23.88 


ND 


78.66 ± 6.60 


15.02 ±2.10 


2227.80 ±145.12 


ND 




5.26 ±0.1 5 


W23 


11 3.57 ±6.72 


ND 


5.57 ±0.85 


34.84 ± 3.75 


73.67 ± 10.12 


ND 




2.56 ± 0.66 


^vearage of three determinations. 
2 Not Detected. 
















within a group and dissimilarities between different groups. 
CA was performed on the data using Ward method and 
squared Euclidean distance. A dendrogram was produced 
by cluster analysis, shown in Figure 2. Seven studied heavy 
metals were classified into five groups based on spatial simi- 
larities and dissimilarities. 


Principal component analysis and factor analysis 

PC A reduces the dimensionality of data by a linear com- 
bination of original data to generate new latent variables 
which are orthogonal and uncorrelated to each other 
[32]. Prior to PCA and FA analysis, the raw data was 
commonly normalized to avoid misclassifications due to 


Table 4 Summery statistics of heavy metal content in water samples (ug/ L) analysis 














Fe 


Co Ni 


Cu 


Zn 


Cd 


Pb 


Detected ( 


%) 




100 


47 78 


100 


100 


65 


100 


Min. of the detected concentration 




75.90 


ND 1 ND 1 


6.59 


27.79 


ND 1 


0.74 


Max. of the detected concentration 




339.75 


99.82 84.15 


65.31 


2227.80 


14.87 


12.45 


Mean of the detected concentration 


168.92 


17.80 24.21 


24.52 


265.12 


2.44 


4.45 


Standard deviation 2 




77.95 


32.33 25.21 


15.50 


456.66 


3.68 


3.14 


MCL (based on WHO) 




300.00 3 


70.00 


1000 


3000 


3.00 


10.00 


Percentage of samples containing metals > WHO (%) 


13.04 


8.70 


0 


0 


17.39 


8.70 


MCL (based on ISIRI) 3 






70.00 


1000 


3000 


3.00 


10.00 


Percentage of samples containing metals > ISIRI (%) 




8.70 


0 


0 


17.39 


8.70 



^ot Detected. 

2 Standard deviation for heavy metal concentration in all samples, 
institute of standards and industrial research of Iran, 1997. 
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Table 5 Spearman correlation coefficient (r) of heavy metals in the sampling stations 





Fe 


Co 


Ni 


Cu 


Zn 


Cd 


Pb 


DO 


Co 


0.019 
















Ni 


0.129 


0.342 














Cu 


0.007 


0.008 


-0.025 












Zn 


0.121 


-0.433 a 


0.593 b 


0.084 










Cd 


0.291 


0.027 


0.679 b 


0.046 


0.551 b 








Pb 


0.283 a 


-0.054 


0.135 


0.583 b 


0.168 


0.374 a 






DO 


-0.079 


0.273 


-0.681 b 


-0.043 


-0.518 a 


-0.436 


-0.250 




Depth 


0.516 a 


0.519 


-0.258 


0.003 


-0.133 


0.155 


0.084 


0.136 


Dist. Ind. 


0.266 


0.027 


0.234 


0.029 


0.080 


0.261 


0.241 


0.063 



Correlation is significant (a) at the 0.05 level and (b) at the 0.01 level. 



the different order of magnitude and range of variation 
of the analytical parameters [30]. The rotation of the 
principal components was executed by the Varimax 
method with Kaiser normalization. 

Four principal components are obtained for heavy 
metals through FA performed on the PC A. This indi- 
cates that four main controlling factors influenced the 
quality of surface water in the study area. Corresponding 
components, variable loadings, and the variances are 
presented in Table 6. Only PCs with eigenvalues greater 
than 1 were considered. PCA of the whole data set 
yielded 4 data sets explaining 88.92% of the total 



variance. First component which explained 32.02% of 
the total variance is correlated with Pb and Cu. The sec- 
ond component is due to Zn and Co. The third compo- 
nent is a location for only Ni metal. The latest extracted 
factor is related to Fe and Cd. 

Discussion 

According to results mentioned in Table 4, all of the 
samples contained Co, Cu and Zn inferior to the values 
specified by related MCLs. In contrast, in 13.0% and 
8.7% of the samples the amounts of Fe and Ni, respect- 
ively, were above WHO MCLs. The amount of cadmium 



Rescaled Distance Cluster Combine 

10 15 20 



15 
_L_ 



-J- 



25 
— L_ 



Pb3 



Group 1 



Cu 4 



Group 2 Zn 1 



Group 3 



Cd 2 



Fe5 



Group 4 Co 7 



Group 5 



N6 



Figure 2 Dendrogram of heavy metal concentrations of water samples. 
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Table 6 Rotated component matrix of four-factor model 3 



Component 





1 


2 


3 


4 


Fe 


0.195 


0.036 


0.267 


0.881 


Co 


-0.129 


-0.867 


0.047 


0.142 


Ni 


0.044 


-0.015 


-0.974 


-0.027 


Cu 


0.960 


0.151 


0.070 


0.144 


Zn 


0.047 


0.858 


0.088 


0.078 


Cd 


0.042 


-0.140 


-0.504 


0.808 


Pb 


0.961 


0.047 


-0.123 


0.078 


Eigen value 


2.241 


1.708 


1.216 


1.059 


% of total variance 


32.020 


24.405 


1 7.366 


15.132 


% Cumulative of variance 


32.020 


56.425 


73.791 


88.924 



Extraction method: Principal component analysis. Rotation method: Varimax 
with Kaiser Normalization. Rotation converged in 5 iterations. 



found in 34.8% of the samples was lower than the detec- 
tion limit of the DPP method, but 17.4% contained the 
metal ion superior than the ISIRI and WHO MCL. This 
is of concern because cadmium has carcinogenic proper- 
ties as well as a long biological half life leading to 
chronic effects as a result of accumulation in liver and 
renal cortex. It can also cause kidney damage as well as 
producing acute health effects resulting from over ex- 
posure to high concentrations [20]. 

Due to possible long term effects of chronic exposure, 
the presence of lead in drinking water is crucially im- 
portant for public concern. Although all of the samples 
included this metal, 8.7% of them contained lead ions 
above the levels proposed by WHO and ISIRI MCL. 
Overall average concentration of heavy metals in water 
samples varies as Zn > Fe > Cu ~ Ni > Co > Pb > Cd. The 
results reveal that the amount of heavy metals depends 
on the sampling locations. 

As shown in Table 5, a close relationship between the 
couples Fe/Pb, Cu/Pb, Cd/Pb, Co/Zn, Ni/Zn, Ni/Cd and 
Zn/Cd states a probable common source of the couples. 
A further statistical investigation was performed by 




Figure 3 Component plot in rotated space for heavy metals (Factor loadings, factor 1 vs. factor 2 vs. factor 3, Rotation: varimax 
normalized, extraction: principal components). 

J 



Zamani et al. Iranian Journal of Environmental Health Sciences & Engineering 2012, 9:29 
http://www.ijehse.eom/content/9/l/29 



Page 9 of 10 



testing the correlation between the determined concen- 
tration of heavy metals and the distance of the sampling 
site from NILZ Company. The calculated correlation 
(Cu/Pb and Zn/Cd) can confirm the significant effect of 
NILZ Company activities as a main source of heavy 
metal contamination observed in the investigated 
groundwater samples. In addition, close correlation be- 
tween Fe and depth of the wells (0.52) suggests that this 
metal is totally of pedo-geochemical source leached from 
the upper soil layers. 

Cluster analysis allows identification of five clusters or 
groups of associated metals (Figure 2). On the basis of 
similarities found for group 1 (Pb, Cu), one can suggest 
the anthropogenic origin of the contamination sources. 
The presence of iron in group 3 (Cd, Fe) notifies, prob- 
ably, mixed anthropogenic and pedo-geochemical source 
of the metals presented in this group. Therefore Zn, Co 
and Ni were located in single member groups. 

Also according to Table 6, Component 1 is attributed 
to lead and copper with positive sign. These elements 
are important byproducts of lead industries indicating its 
anthropogenic sources. Component 2 reveals 24.4% of 
the total variances are positively loaded with Zn and 
negatively loaded with Co. Component 3 shows that 
17.4% of the total variance is positively loaded with Ni 
and it can be represented by oil industries activities near 
the NILZ Company. Component 4 explains 15.1% of the 
total variance, is positively loaded with Cd and Fe. 

The heavy metal grouping has been explored in the 
plot of the first three principal components generated 
from these parameters (Figure 3). The low correlation 
found for the studied metal ions in the four components 
defined by FA, suggests both anthropogenic and pedo- 
geochemical sources for the metal contaminations. 

Conclusion 

Overally, the present study has shown that the ground- 
water source within radius of 5 km from National Iran- 
ian Lead and Zinc Company (NILZ) in Bonab Industrial 
Estate (Zanjan province-Iran) is contaminated by iron, 
cobalt, nickel, copper, zinc, cadmium and lead. This can 
be considered as a menace for people who daily intake 
the corresponding waters, planted vegetables and food 
crops irrigated by the same water source. The higher 
health risk comes from those elements which are present 
at higher levels than announced by WHO and ISIRI not- 
ably lead, nickel and cadmium. Multivariate statistical 
techniques have shown correlations and similarities 
among the investigated heavy metals and classification 
of these ion groups. Cluster analysis has identified five 
clusters among the heavy metals. The statistical investi- 
gations reveal the pollution sources influencing water 
quality in the study area as anthropogenic (with a very 
high contribution of NILZ Company) and pedo- 



geochemical for Fe, Cu. The results suggest a significant 
risk to the population of Zanjan city and its neighbor- 
hoods given the toxicity of the studied metals and the 
fact that this aquifer by far is the main source of their 
drinking water and irrigation. This study has also high- 
lighted the need for further research and regular moni- 
toring, in order to determine the permitted levels of 
metals in the studied aquifer. 
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